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I. INTRODUCTION 


A. WHY SWITCHED CAPACITORS? 

Modern active filter design 1s based mainly on RC filters. Practical RC filter 
realizations require in most cases large-valued capacitors, and high precision resistors 
and capacitors components to achieve accurate time constant. These two features 
created a major obstacle for IC implementation of RC filters. 

The MOS integrated-circuit technology had found wide usage in industry because 
of its superior logic density, as compared to that achievable with bipolar technologies. 
A unique property which is held by the MOS integrated-circuits is the capability to 
store charge on a node for a short period of time, order of milliseconds, and to sense 
this stored charge continuously and nondestructively. This property was first used in 
dynanuc random-access memories and dynanuc logic. In the last few years the analog 
signal processing development gained a lot of momentum using the property mentioned 
above. 

As it will be explained in the next section of this chapter, the realization of the 
RC time constant can be well controlled in MOS technology by determining clock 
period and the capacitor ratio. This idea made it possible for active filters to be 
implemented in IC form. Furthermore, MOS capacitors are nearly ideal, with very low 
dissipation factors and good temperature stability. As an economical approach, 
sampled-data filters can be fabricated using memory-like NMOS and CMOS 
processing. This analog and digital circuitry can be placed on the same chip. For the 


reasons presented above switched capacitor filters became attractive and useful. 


B. NOTATION 

In general sense, a switched capacitor network (SC) would be made up of 
switches, driven by an external clock, capacitors, and OP-amp if it is needed. 

It is appropriate to identify some of the notations and symbols that will be used 
throughout this thesis. The switch symbol will always be shown open. It is said that 
the system has an n-phase clock, if the period of clock, T, has n segments 
corresponding to each phase. For simplicity, a two-phase clock is going to be used. 


Each switch associated with the proper clock phase will be designated by the symbol 


®, 


Where k is an integer representing the number of clock phases in the system. Typically 
a switch is closed only once during the clock period. It is very important in SC 
networks to ensure that at any time, no two switches driven by different clock phases 
will be closed simultaneously (short circuit). Therefore, each phase of the clock has to 
have a less than 50% duty cycle. This is known as the nonoverlapping clock. Bipolar 
two phase nonoverlapping clock is shown in Figure 1.1. A convenient complementary 
notation for the clock phases is denoted by odd defined as ®, and even, defined as ®,. 
Thus, the sampled data waveforms can be stated as the sum of their odd and even 


components, e or o. 





Figure 1.1 Waveforms of a Two-Phase Clock. 


Ei 


If we want to express the transfer function of a SC network, we can write at least 
four transfer functions corresponding to even and odd phase. A useful notation for the 


transfer function 1s shown in Equation 1.1; where 1 andj can be either e or o. 


Vq/(Z) 


v.12) (eqn 1.1) 


H'(z) = 


Since the switched capacitor networks have sampled data character, they can be 
treated similar to digital filters and analized in z-domain. However, SC filters are 
analog networks; thus, the analog concepts of impedance and loading, which are 
absent in digital filters, are retained. To solve this problem and to be able to apply 
network theory, z-domain equivalent circuits, are used, (also known as building 
blocks). 


C. SWITCHED CAPACITOR EQUIVALENT RESISTORS 

In this section realization of continuous resistor will be denoted in several 
methods using switched capacitor networks. The main goal here is to have the ability 
to replace the continuous resistor using an equivalent SC realization. Several methods 
are developed for continuous resistor realization. Four of these realizations are 
discussed below. First, is the parallel switched capacitor resistor realization. In this 


method, the circuit in Figure 1.2 is considered. 





Figure 1.2 Phe Parallel SC@ Resisterealizanorne 
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Initially, both switches are open and the capacitor is completely discharged. In 
the odd phase, M, is closed and capacitor will be charged to V, this can be expressed 
in the following 


Q(tg + 1/2) = CV, (eq ral. 2) 
Because of the nonoverlapping clock, first @, will be open and then even phase 


will start, D, will be closed, and C will be charged to V,. In this phase, total charge on 


the capacitor is 
One 1) = CV, — CV, (eqn 1.3) 
We note that the final charge on the capacitor during one phase period is not 


necessarily equal to the charge flowing past the voltage sources during that period. 


The current can be written as 


i= dq dt (ean 1.4) 
dq = CdV (eqn 1.5) 
eels Vv, } (eqn 1.6) 


on the average 
_ eo5— Vv) 


(eqn 1.7) 
T 





= we (eqn 1.8) 
ie 


This vields the following relationship for parallel realization 


R = T/C (edna) 


where T is the clock period. 

Second, let us consider the circuit shown in Figure 1.3 called the series switched 
capacitor realization of the continuous resistor. Using the same previous calculations 
for parallel SC resistor realization we can obtain Equation 1.9 as equivalent resistance 


for this type of realization. 





Figure 1.3 The Series SC Resistor Realization. 


A third realization is a combination of the parallel and series realization. This 
type of circuit is shown in Figure 1.4. Initially, again, both capacitors are discharged 


and switches are open. In the odd phase, ®, is closed and C, is charged to V). 


ONG, > 2) ee (eqn 1.10) 
In the even phase, ®, is open and ®, is closed respectively. 

Q,(tg+T) = C(V,— Vy = GW, -¥) (eqn 1.11) 
consequently 

Qigjt l= CVs) (eqn 1.12) 


Equation 1.12 in fact is the component of Q,(tp+T) with the minus sign. If we 
continue to take the other phase, which is odd phase, we can write the following 


equation 


Qi (tg + 3Ti2) = CV, — V5) (eq le 13) 


So far, we completed one clock period, from T/2 to 3T/2. Thus, we calculate the 
total charge going into voltage sources. 
— ee a ei Ww; 
Q) = ypPydt = Paydet fide 
eee aie e y(t, + 3; 2) 
where u= tyt T/2, v= tpt I, w= tyt 3T/2. On the average 


Cee Civ Ve 
L = =< (eqn 1.14) 


y= (eqn 1.15) 








Figure 1.4 The Combination of the Parallel and Series Realization. 


This realization is known as series-parallel realization of the continuous resistor. 
As it was seen from the Equation 1.15 the equivalent resistor depends on C, and C,,. If 


the capacitor values are chosen equally, then Equation 1.15 becomes 
R = T/2C (eqn 1.16) 


Since twice the charge is transfered from one end to another, the resistance value 


is reduced by 2, as compared to the previous realizations. 


There is another realization of the continuous resistor which has indeed the same 
result as the series-parallel realization. This configuration 1s shown in Figure 1.5 and 


known as the bilinear switched capacitance realization of a continuous resistor. 





Figure 1.5 The Bilinear SC Realization of a Resistor. 


In this approach, all switches change position twice every clock period. 
Therefore, actual clock period applied to circuit is T/2, rather than T. We can put T/2 


in place of T in Equation 1.16 and obtain 
R= T/4C (eqn 171) 


In order to show some practical differences of these realizations, some outcomes 


of the experiments will be discussed and displayed in the next sections. 


D. CHARGE CONSERVATION ANALYSIS 

This method is a slightly different tvpe of application of Kirchoffs current law, 
Where charge q is used rather than current. [In general, because of the two phase 
clocking operation, a pair of charge equations can be obtained, which characterize the 
charge conservation condition at a particular node for all sample instants. A pair of 
charge equations consist of both equations for the even sampling instants and for the 
odd sampling instants. These are written as follows 


for even clock phase (®,) 


qc) Sq ae eat 
for odd clock phase (®, ) 
q”(t) = “Gee es >t 
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Where q,(t’) is the charge held at one particular node at reference time t'’, dit) is the 
charge at the particular node from the previous phase period, called memory charge, 
and q(t) is the charge injected at that particular node, called contribution charge. 

SC networks are characterized mainly in terms of charge-transfer’ functions on 
which the discrete-time voltage v(kT) and discrete-time charge differentation Aq(kT) 
are used as variables. 

Assuming a SC network with i capacitors, and N denotes the total number of 
capacitors connected to a node P during even and odd clock phases, charge variations 


can be written as 

AGC) — crip ea veo ee co 1) (eqn 1.18) 
for even n, and i= L250 Nop: 

Ag® (nT) = ),0°,(0T) — 0° (a - 17) (eqn 1.19) 


fo~woed n, and i=1,2,...N_. 
op 


In the z-domain, the above charge equations can be described as follows 


AQ(Z)= LX (2) —Z FAL V2) (eqn 1.20) 
where 1= Lys Nop 
“12 
AQ’ (2) > > jQ,i(2) ~ 2 >, 2, Z) (eqn 1.21) 
where i= 1,....No,; 


Equation 1.18 through Equation 1.21 are going to be used throughout this research. 
E. EXPERIMENTAL DEMONSTRATION OF SWITCHED CAPACITOR AS 
EQUIVALENT RESISTOR 
In this section, the mathematical and theoretical relationship between SC 
networks, of which the output is a sampled analog signal, and the corresponding 
equivalent resistance value found in the previous section, will be verified experimentaly. 
As discussed previously, a two phase clock will be needed for all SC network 


realizations. These two phases can be practically generated from a single clock using a 


simple digital circuit shown in Figure 1.6 This circuit is capable of producing two 
nonoverlapping output clocks that can be used as the two phase clock needed. Output 
of the NOR gate will be high whenever both inputs are low. This ensures that one 
phase will be high when the other is low. To provide a larger margin between the two 
phases of the output clocks (longer nonoverlapping time), two inverters were utilized in 
each feedback loop to create a longer delay. In the lab, CD4049AE CMOS bipolar 
gates were used as inverter, which has approximately 15 ns delay time, and CD4001AE 
CMOS bipolar gates were used as NOR gate, which has 20 ns delay time. Master 
clock, which a is bipolar square wave, was provided by EXACT MODEL 120 
WAVEFORM GENERATOR. The output of the clock circuit is shown in Figures 1.7, 


1.8, 1.9. The nonoverlapping time was found to be 50 ns as expected. 


L MASTER CLOCK 


| CDAOOIAE 


CD4049AE CDS049AE 
CD4049AE 


AOAIAE 
CD CD4049AE CD4001AE 





Figure Ge The lworphasci@lec@incuit 


Two examples of SC resistor realizations were considered. The first 1s the parallel 
switched capacitor realization of a continuous resistance. A simple voltage divider 
network was constructed using this realization to demonstrate how close the 
experimental results are to the computed values. This circuit is shown in Figure 1.10. 
The switches used in this circuit are the CMOS Bilateral Switches CD4066B. These 
devices are most suitable for switched capacitor applications due to their lower channel 


resistance (80 ohms maximum), and will be used in all the experiments of this research. 
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Figure 1.9 Falling Edge. 





Figure 1.10 Experimental Circuit Using Parallel SC Realization. 


The results of this experiment using 50.96 KHz and 101.19 KHz clock frequencies are 


shown in Figure 1.11 and Figure 1.12. 


PARALLEL SWITCHED CAPACITOR REALIZATION OF RESISTANCE 


10° 10° 


RESISTANCE VALUE 
104 


10? 





CAPACITANCE VALUE 


Figure 1.11 Experimental Result of Parallel SC Realization. 


PARALLEL SWITCHED CAPACITOR REALIZATION OF RESISTANCE 


CLOCK FREQUENCY = 101.19 KHZ 


10° 10° 


10° 


RESISTANCE VALUE 





Figure 1.12 Experimental Result of Parallel SC Realization. 


a 


The second example ts the bilinear switched capacitor realization of a continuous 
resistance. This experiment was conducted to compare the experimental results 
obtained using this topology and the theoratical value given by Equation 1.17. The 
results were also used to show the advantage of this topology over the parallel 
realization technique. The circuit used in this experiment is shown in Figure 1.13. 
The result using two different clock frequencies, 51.1 KHz and 101.82 KHz are 
displaved in Figure 1.14 and Figure 1.15. 

In all of the above results, one can notice some deviation between experimental 
and computed values. These are caused by nonideal properties of the switched 
capacitor which are going to be discussed later in this section [Ref. 1]. 

These experimental results show that different SC realizations can be used to 
simulate resistors when switched capacitor value is within certain range depending on 
the clock frequency. Otherwise the use of equivalent resistance can lead to gross errors 
[Ref 2): 

Altough the bilinear realization yields a better approximation, there is an 
important disadvantage which is known as parasitic capacitance. This idea may be 
shown in Figure 1.16. The parasitic capacitance, Cg, is the capacitance which 1s 
virtually placed between one or both plates of the switched capacitor and ground. The 
value of this parasitic capacitance can reach 10 percent of the switched capacitor 
ecien3) 

Because of this parasitic capasitance, final equivalent resistance value were found 
to be lower than expected as seen appearently in Figure 1.14 and 1.15. For this reason 
the value of C can be selected a few percent lower than the nominal value to neglect 
some of effects of Cg. When C is lower than the computed value, according to 
Equation 1.17, R becomes larger and matches the computed line in Figure 1.14 and 
Figure tis 


F. NONIDEAL PROPERTIES OF THE SWITCHED CAPACITOR 
The use of CMOS FET as switches introduces some secondary effects to the 
switched capacitor network. These effects may be classified as follows [Ref. 1]. 
1. Clock feedthrough. 
2. Offset error and noise. 
3. Nonlinear P-N junction capacitance. 
4. Incomplete transfer of charge. 


These secondary effects will be studied throughout this research. 


Re 
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Figure 1.13. Experimental Circuit Using Bilinear SC Realization. 
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Figure 1.14 Experimental Result of Bilinear SC Realization. 
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Figure 1.16 Appearing of the Parasitic Capacitance. 
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Il. THE BILINEAR TRANSFORMATION AS A SC FILTER DESIGN 
TE OE 
A. GENERAL 

Since the switched capacitor filters are sampled version of analog filters, the main 
goal of this chapter is to set up a bridge between continuous and discrete time 
domains, and use this as a design tool. 

Since filters are typically specified by frequency-domain requirements, it is 
convenient to have a mathematical expression that allows to transform rational s- 
domain transfer functions to rational z-domain transfer functions. Such an expression 
should possess two qualities: 

1. Stable s-domain transfer functions map into stable z-domain transfer functions. 
2. The imaginary jw-axis of the s-plane map onto the unit circle of z-plane. 
Item 1 ensures that transformed z-domain transfer functions will be stable. Item 2 
ensures that not only will they be stable but that the shape of the gain response can be 
preserved. 

From digital signal processing environment, there are several major mapping 
procedures which are used to transfer the information from s-domain into the z-domain 
or vice versa. They are as follows 

1. Backward Difference (BD) Transformation 

2. Forward Difference (FD) Transformation 

3. Lossless Discrete Integrator (LDI) Transformation 
4. Bilinear Transformation 

By setting s= — @ + jw and evaluating z, it can be determined whether items | 
and 2 are satisfied. The backward difference transformation maps the jw-axis inside 
the unit circle. Thus, item 1 is satisfied; however, item 2 is not. This mapping is stated 


as follows 
lela SL) (eqn 2.1) 


Appliying the same procedure for the forward difference transformation, it can be 
showed that the FD maps the jw-axis into a straight line outside the unit circle. Thus. 
the forward difference transformation satisfied neither item 1 nor 2. This 


transformation can be expressed as follows 


2D 


Ze sil (Comeze2) 


It can be concluded that the forward difference transformation maps s-plane poles and 
zeros with high quality factor into z-plane poles and zeros with higher quality factor. 
The resulting z-plane poles may be unstable. One can determine that the FD is 
unsuitable and that the BD is marginally suitable (for wI < < 1) for transformation. 
The lossless discrete integrator (LDI) transformation and the bilinear 
transformation satisfy two qualities.But, the LDI transformed transfer functions being 
realized require twice as much hardware as that required to realize bilinear transformed 
transfer functions. In the previous chapter, the results of the experiments showed that 
the bilinear realization of the continuous resistance would appear better than others. 
The bilinear realization is also very related to the bilinear mapping procedure that will 
be explained later in this chapter (Section E). For these reasons, the bilinear 
transformation will be deeply studied and vital points underlined so that one can use 


these points as a design aids. 


B. BILINEAR TRANSFORMATION 

The bilinear transformation 1s a method which is based on numerical analysis 
techniques. It solves the difference equation, which represents the system, according to 
the trapezoidal rule. In order to illustrate the technique let us follow an example. As a 
notation, x(t) is used for input signal of an analog system and y_(t) as the output 
signal of the system. Let us assume that these two signals are related to each other by 


the first-order differential equation 


dy, (t) 


b 
1 dt 





+ Boy (t) = agx(t) (eqn 2.3) 


The corresponding analog system transfer function 1s 





a 
H(s) = 0 eqn 2.4 
(s) b s+ by (eq ) 


y,(t) can be expressed as follows 


Salt ma 


> 


Y,(tg) (eqn 2.5) 


In this step, if the signal is sampled every T second and if t=nT, th > (n-1)T, then 
y(nT) = tyr) lya(r)de + y ((n-1)T) (eqn 2.6) 
If the integral process is evaluated by the trapezoidal rule we obtain 
y4(nT) = y,((n- 1)T)+ T/Jy’,(nT)+y/,((n— HT) (eqn 2.7) 


From the original differential equation 


bo Ay 
a) = eee 
substituting into Equation 2.7 we obtain 
pee ao) = 12 -baib ay (nt) ylides) 0) a/b, (x (mn) x (n-1)T))). 


Taking the z-transform and solving for H(z) gives 


eo 2.8 
De a or Se 


am itz! 8 





If H(z) is compared with H(s), relationship between s and z can be obtained easily and 


Written 


H(z) = H(s) | oy {eqn 2.9) 


This can be shown to hold in general since an N"-order differential equation can be 


written as a set of N first-order equations of the form of Equation 2.3. 


2 1-z'! 
= => (eqn 2.10) 
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Or 


1+(1/2)s 
Z = ——_—— (ciieeduly 
F=( 1; 2s 
That is, the invertible transformation of Equation 2.10 is recognized as a bilinear 


transformation. Setting s=a@+y]w give us 


— 1+(T/2)(a+ jw) 


= e aly 
1—(T/2\a+jw) oe 2 


by letting z = rei? we obtain 


J/((2/T) +4)? +w? 


r= See and 
J (21) = a+ w? 


ay WwW 
Se an eqn 2.13 
Ita? ae ee 


5 Ww 
0 = tan 


From Equation 2.13 the bilinear transformation maps 


1. The right-half of the _s-plane onto the region exterior to the unit circle where 
zZ|= 1 of the z-plane. This comes from the fact that @>0, r>1. 


2. The imaginary axis of the s-plane onto the unit, circle, |zZ|=1, of the z-plane. 
This comes from the fact that @=0, z= exp|j2tan (wT/2 


3. The left-half of the s-plane onto the interior of the unit circle, |z|=1, of the z- 
plane, because for a<0, r> 1. 


From item 2 several other important observations can be made. For «=0 and w=0, 
z=1e which indicate that the s-plane origin maps to the point (1,0) in the z-plane. 
For a=0 and w>9 SX, z= le® and for @=0 and w> 0O,zZ= te) which indicate that 
the positive and negative w-axis in the s-plane map onto the upper and lower semi- 
circles where |z|=1 in the z-plane. Obviously, a stable analog filter produces a stable 
digital filter. 

Since the entire imaginary axis (s=jw) in the s-plane, maps onto the unit circle 
(/zZ}=1) in the z-plane, nonlinear relationship between analog and digital frequencies 
can be easily seen. The distortion introduced here is known as the warping effect and it 


is illustrated using Equation 2.10 for s=jw and z= if, 
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ey 
W=--= 
meee 


e 


2 iD: 2(¢i8 2 ae ot =) 





2 sin(®.2) 


T © cos(8/2) 


(eqn 2.14) 
and thus the imaginary axis in the s-plane is related to the unit circle in the z-plane 


through the relation 
me.) tan(0) 2) (eqn 2.15) 


where 9 is the digital frequency and w is the analog frequency. 
To obtain H(z) from H{(s) is rather tedious and time consuming. This 
transformation thus is generally done by the computer. For a second-order transfer 


function a computer program has been written and presented in Appendix A. 


C. THE EFFECTS OF HIGH SAMPLING RATE 

In the mapping process, generally the s-domain specifications of the filter are 
known. But, there might be some ambiguity if it is the case of obtaining a function 
H(s) from a given function H(z). However, when the sampling rate is much greater 
than the frequencies of interest, various assumptions can be made to obtain 
approximate relationships between the two domains. 


For small |s|T, Equation 2.11 can be written 


(sT)? (sT) 
z= 1+s le oo (eqn 2.16) 


For high sampling rates, the second-order effects can be ignored, so Equation 2.11 can 


be approximated as 


Je eS | (eqn 2.17) 


Equation 2.17 shows that as the sampling period T is decreased, or sampling frequency 
increased, z approaches the point z= 1. 


The characteristic equation of the second-order analog filter can be expressed as 





W 
s7+ 0 : st Wee = 0 where Q is the quality factor (eqn 2.18) 


The solution of the characteristic equation, that 1s the poles of the system can be 


described as 


W \ a 
Ss . = — —14)—2 /407-1 eqn 2.19 
12 20 Sov Q (eq ) 





On the other hand, the characteristic equation of the second-order digital filter can be 


written as follows 
27> — 2rCosO0z+1r*? = 0 (eqn 2.20) 


using the approximation obtained in Equation 2.17 


21> I+s) 41 consequently 
Se re #18 
_ iW Tie 
r= J/1-(W,T/Q)+(W,T) (eqn 2.21) 
W_T/4Q?= 1 
8 = tan”! Bala we | (eqn 2.22) 
Ome Wied 
if Equation 2.21 is solved for Q 
WT 
Q 4 (eqnigzsz 3) 


~ (ln +n+(wo Ty 
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if Q is greater than unity as well as assuming high sampling rate, Equation 2.23 can be 


approximated as 


W.T 


—_-2___. 2.24 
21 —1)+(W TY ue 


~w 


In the case of high sampling rate warping effect can be neglected and §@~W_ T. Thus 


Equation 2.24 becomes 


8 
20 =r) 





Q~ ore (rea esr (VV ally) (eqn 2.25) 


from Equation 2.25, r can be approximated as 
r=1-— (eqn 2.26) 


and Cos@ can be written as 


G2 
Cos8 = 1- ai (eqn 2.27) 


Note that Equation 2.21 through Equation 2.27 are valid at high sampling rates and 
for large value of Q. For simplicity, if the characteristic equation of the z-domain 
transfer function is considered as 

2 ae 

Tae Ze. = -0 


eee, lb) b, 


Wo) = Kk 
- Ale) 1D 


(eqn 2.28) 


For bilinear transformation, 


7 ele by bs )(1 = 6) 65) 
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Ze >) vas 


oi 


where K is a constant involving a frequency prewarping [Ref. 4] and is given by 


K is ( 70 
= ——___ eqn 2. 
tan(@/2) : 


When high sampling rate condition is satisfied, 0/2 < < 1, K simply becomes 


K = 2/T (eqn 2.31) 


D. ELEMENT TRANSFORMATIONS 

In order to treat SC filters in z-domain, some concepts which don’t exist in 
digital filters should be considered. These concepts are basically associated with 
passive circuit elements, which are inductances, capacitances and resistances. 

If the bilinear transformation is used directly to obtain the admittances, one can 
use the results as an equivalent elements which are derived as follows. This method is 
usually known as Building Block Approach. 

1. Capacitor 

If the relationship between current and voltage across a capacitor is written as 


(t= ee is, dq(t)_ _ dv(t) 


dt 
by taking Laplace transform of both sides 


sQ(s) — q(0) = cf{sV(s) — V(0)] where, L[q(t)] = Q(s) and YH[v(t)] = V(s). 
q(0) = v(O) = 0 as initial conditions. 


Q(s)/V(s) = C 


using bilinear transformation 


Y=—-==C (eqmer: 32) 


Equation 2.32 can be considered as admitance, which is not the way it 1s 
described. This new definition comes from the fact that the analysis of the SC circuit 1s 


done by using charge variation equations. 
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2. Resistor 
If the same procedure followed for capacitor is applied for resistor, the 


following equations can be easily written. 


Since it) = Gv(t) Where G = I/R 
da(t 
dt 


taking Laplace transform of both sides 
sQ(s) — q(0) = GV(s) where q(0) = 0 
Qs) _ G 


V(s) S 


using bilinear transform vields 


Yop = — = — = ee 
R Siz} S 2 1-z'! sean 


3. Inductor 


The instantaneous voltage and current for an inductor L are related by 





di(t) 
v(t) = L— 
(t) rs 
d-q(t 
v(t) = L a2 
taking Laplace transform 
V(s) = L{s*Q(s) — sq(0) — q'(0)] where q(0) = q‘(0) = 0 
O(s) | 7 a 
— = -— > using bilinear transformation yields 
(s) Ls 
Ce) UG 


iS epee ee eqn 2.34 
L vz) Ls? AL (1-27!) 3 


To express the single floating capacitor shown in Figure 2.1, charge variation 
equations, Equation 1.20 and Equation 1.21, can be used. 

AQ®,(z) = CVE,(z)~ Cz" V9 (z) ~ CVE, (2) + Cz" 2 V°(z) 

Oe (Ze CV (Zz) — CZ AVE (2) — CV 32) ae BN) 

AQ®,(z) = CV°,(z) — Cz" #V°,(z) ~ CV®,(z) + Cz" 7-V° (2) 





Figure 2.1 Floating Capacitor. 


AQ®,(z) = CV°,(z)— Cz" #V°,(z) - CV®, (2) + Cz" V8 (2) 


From these equations 


AQ*(z)+ AQ? (z) = C(L— 2" )[V* (2) + V° (2)] (eqn 2.35) 
— C(1—2"?)IV®,(z) + V°,(z)) 


and 


AQ®,(z) + AQ°,(z) = C(L-z" )[Ve (2) + v°.(2)] (eqn 2.36) 
= C2122) it ven 


Since AQ(z) = AQ%(z)+AQ%z) and 
V(Z) "Se Vista) een) 


Equation 2.35 and Equation 2.36 can be written as 


AQ (z) = C(l=z 7)V,(2)-C =z 2V2) (eqn 2.37) 
AQ.(z) = C(1-2°7)V,(2)- C1 - 2° 7) V (2) (eqn 2.38) 
From Equation 2.37 and 2.38, admitance Y, appears as Y, = C(1—~z /2) 


° - as ip ate cea / a a 
Since z! = eW¥! then 2? = ei¥l2 even though the period of the sampling 
decreases to half of the orginal value. Then 


; 
Zs z”2 can be written. 


we C(l=@)) (eqn 2.39) 


Since Equation 2.39 is realizable with simple SC networks, we can compare this, ie 
with Equation 2.32 and consider (1—(z’)') factor as a scaling factor [Ref. 5]. 
Therefore Equation 2.32, Equation 2.33, and Equation 2.34 can be written as 


= oe el 
= (lz) 





(eqn 2.40) 
Yp = GT/2(1+2z"') (eqn 2.41) 
ie (letza!)2 cis, 
= — CCiiez: 
EL 4p d@—2)) " 


This would be a good place to revisit the bilinear switched capacitance realization of 
the continuous resistor shown in Figure 1.5. It is also known as the floating bilinear 
resistor. Charge equations in z-domain are 

AQ®(z) = CV*(z)+Cz V9, (z) — CV*,(z)— Cz"? V°,(z) 

AQ° (z) = CV° (2)+Cz V2, (z)— CV°,(z)— Cz? V°,(z) 

RAZ) — CV(z)+ Gz) 7V°(2)- GV" (z)- Cz AV" (2) 

AQ°,(z) = CV°,(z)+ Cz" “*V*,(z) — CV? (z)— Cz" V8, (2) 


AQ (z)t AQ, = C(1+z°)V,(z)— CU +z 7)V,(z) (eqn 2.43) 


AQ®,(z) + AQ®,(z)= C(1 +2 7)V,(z)— CU+z 7)V, (2) (eqn 2.44) 
Equation 2.43 and Equation 2.44 yields 

Yp = C(l+z'”) (eqn 2.45) 
Since (z’)! = z 12 this can be written as 
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ape os C1 az ae (eqn 2.46) 


if Equation 2.46 is compared to Equation 2.41 
Yp = (GT/21+(zy) = C+ (zy) 


G = 2C/T 
T’ 
R= — andT = 1/2 
2C 
R 2 ( 2.47) 
=— eqn 2. 
aC dq 


which is the same as Equation 1.17, as expected. 


E. INTEGRATORS USING BILINEAR TRANSFORMATION 

Since the integrators are widely used circuits in analog signal processing, their 
transfer functions will be investigated and compared with the results obtained in the 
lab to decide how much error would be introduced using bilinear transformation. 

1. Lossless Integrator 

An analog lossless integrator is shown in Figure 2.2. If the resistor is replaced 
by the bilinear SC resistor equivalent circuits presented in Chapter I, Figure 2.3 is 
obtained and obviuosly 
r 


C 
dR, 


Figure 2.2. An Analog Integrator. 





Figure 2.3 Bilinear SC Integrator. 


if charge variation equations are written in the z-domain, 
Gave ©, Vz) CN cz me) 
can be obtained for odd phase and 
e te 0 4 o -Vy e 
C,V a eon. 3(2)z mee 220 oe ee 2) 


can be obtained for even phase. Knowing that 
me) = Vv 2+v@ 


beige V2) V2) 
from Equation 2.48 and Equation 2.49 





Hg) = 22 2 
V (2) C27) 


LY, oa 
Since z’? = (z’)7! the transfer function H(z) is rewritten as 


ee zor 
oe 1-77) 


Transfer function for a simple analog integrator is 





l ali 
iS) = where R, = —— , asin Equation 1.16 
R,C,S ue 


(eqn 2.48) 


(eqn 2.49) 


(eqn 2.50) 


(ogi 2.51) 


aE; 


Cans 


H(s) = - 





(eqnizs2) 


If Equation 2.51 is compared to Equation 2.52 
2 = 
S$ =  o_ 
ee ' , 
which is nothing more than bilinear mapping. 





The lossless integrator would not be implemented due to the fact that the 
Operational amplifier requires some dc feedback. The computer simulations of both 
H(s) and H(z) are shown in Figure 2.4 and Figure 2.5. These figures demonstrate 
clearly the nonlinear relationship between s-domain and z-domain due to the warping 
property introduced by the bilinear transformation. 

2. Lossy Integrator and Experimental Results 

As mentioned earlier, the lossless integrator can not be realized in practice. 
Instead, the lossy integrator is considered in practical SC circuit. The lossy bilinear 
integrator is shown in Figure 2.6. The analog lossy integrator’s transfer function H(s) 


can be written as 


V4(s) ae ] 


H(s) = > ——-—_____ 
(s) R,C.(s+ 1(R,C,)) 


(eqn 255) 


where R, = T/(4C,) and R, = T/(4C;). 
If a high sampling rate is used the warping effect can be ignored. The bilinear 


transformation yields the following 


C (l+(zy4) 


H(z’) = = 2 
C, 2™t(zy}+CyC,(1- (zy 


(eqn 2.54) 


As an experimental set of values in the lab, 


Com elon 
Carre 
C, = 4.2 pf 


were choosen and CD4066 CMOS switches and LM74I1 Operational Amplifiers were 
used. H(z’) was obtained by the computer program given in Appendix A, using the 


following equation. 


38 


Re ee oe a la, 
FREQUENCY RESPONSE MAGNITUDE 


MAGN T TUDE 


1000.0 180.0 2000.0 2500.0 3000.0 3800.0 
FREQUENCY Hz 





Figure 2.4 Bode Plot for an Analog Integrator. 
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Figure 2.5 Magnitude Response for the Bilinear SC Integrator. 
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Figure 2.6 The Lossy Bilinear Integrator. 


H(s) = — 1.066 10° (eqn 2.55) 


s+ 1.9841 102 


| 0.5195(z’)* + 1.0392’ + 0.5195 
H(z’) = -— —————————_ with f= 50.8 KHz (eqn 2.56) 
(z')* + 0.0193z’ — 0.98066 : 


The computer simulation results were displayed in Figure 2.7, Figure 2.8, Figure 2.9, 
and Figure 2.10. Using building block approach, the transfer function of Figure 2.6 can 


be readily written as follows 


| Cu acm) 
H(z) = + (eqn 2.57) 
Cy (la Sy ©) a ae ez 
Using the previuos capacitor values H(z’) becomes 
— 0.5243z' — 0.5243 
H(z’) (eqn 2.58) 


1.009742’ — 0.99026 


Equation 25S 1s shown im Figure Zen leandeirounes2 ae 

Finally Figure 2.13 shows the output of the practical circuit of Figure 2.6. The 
experimental result perfectly matches all the computer simulations outputs of Figure 
2.7 through Figure 2.12. Further more, a smoothing low-pass filter is connected to the 
output of the practical circuit to get rid of the clock feedthrough problem mentioned in 
Chapter fe 


40 









R 1 









Sel ae 451 pF ome dil 
FREQUENCY RESPONSE MAGNITUDE 







cece ce Oee ce yoo eee eee ee ss SOS OS OHH OS FOS COE OS COED EOOHODD SOOO ES OO DESOTO TESS ED EEE OSE OOE HSER eoeeseeee ee yeerseeaeenD 


eeeece cocobecouccscosobesosvccsvashescesscssscbscseccsscsstsccocsssossisccccosessolossccscccesiancserecces Pocceccscces 
‘ 1 1 


See Sees eee Fees SF SSHSHHTHESHHEEHSHUFSHSHSSFSHSSHS SHS SSSSEFSHSSSHSHFEFSHSSPFSSSSHSSSSST SHC SSSFHFESFFSSSHEHEHSHSHEHSHFFSHsseesesssvesgegesds 


MAGN ETUDE 


Pee ee ee eo Pee owt EOC O EERE HE OOE TE UOT ECHL OES POE EEE EO DEO OH CODE O ECO PES HE OE OOSE EEE UHOE ES OO EE EOEOCOOOE OODLES ES OE OET OR 





eeeececse ce eee cessseseest sen eeeeeaeates cc cccetecscccccscctonscccccoccsisecccetesscctoccscccsccetocscscssecest seusccscese 


3 1 3 | | a | 3 ; 
0.0 500.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 4500.0 5000.0 
FREQUENCY Hz 






Figure 2.7 Computer Simulation of Lossy Integrator in s-Domain (Magnitude). 
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Figure 2.8 Computer Simulation of Lossy Integrator in s-Domain (Phase). 
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Figure 2.9 Computer Simulation of Lossy Bilinear Integrator in z-Domain (Magnitude). 
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Figure 2.10 Computer Simulation of Lossy Bilinear Integrator in z-Domain (Phase). 
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Figure 2.11 COTS Simulation of Lossv Bilinear Integrator 
Using Building Blocks Approach (Magnitude). 
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Figure 2.12 Computer Simulation of Lossy Bilinear Integrator 
i. ising Building Blocks Approach (Phase). 
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Figure 2.13 Experimental Result of Lossy Bilinear Integrator. 


These results showed that the following design steps can be used to meet the 


given specifications. 
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Definitions of the filter specifications in the frequency domain. 

Construction of the continuous transfer function H(s). 

Transformation from H(s) to H(z’) using bilinear mapping technique. 

Use of building block approach as implimentation. 

Sey construction of switched capacitor realization of the filter. Or after step 


Resistor simulation of the continuous resistance. 
Go to step 5 (taking into account nonideal properties of the SC networks). 


The most important outcome of this experiments is the different approaches 


one can use to simulate continuous resistance using bilinear switched capacitor 


realization. In one approach the bilinear transformation is used through the step I, 2, 


3, 4, 5. The other approach would yield the same results using steps 1, 2, 6, 3, 7. This 


fact is shown in Figure 2.11 and Figure 2.9 respectively. 
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Il. PROGRAMMABLE GIC SC FILTER 


A. GENERAL 

The programmable filter is a very interesting and useful application of the SC 
networks. This is due to the fact that the SC circuit’s frequency response can be 
changed by varying clock frequency. In fact a remarkable feature is that, multiplying 
the clock frequency by a factor @ results in scaling the frequency response (multiplied) 
by the same factor along the frequency axis. 

Most of the active SC networks have limitation over the clock frequency range, 
depending on the operational amplifier’s settling time used in the circuit. In today’s 
technology, \f[OS op. amps. can be designed with 0.5 [sec settling time to reach 0.1% 
of its final value of the step response. Clock rates up to 1 MHz seem feasible. In 
practice, VIOS op. amps. has been designed which settle to within 0.1% of the final 
value in 2 jisec and achieve dc gains greater than 60 dB. Therefore sampling rates of up 
to 250 KHz, can be easily achieved. Considering the Nyquist criteria, almost a 100 
KHz signal can be processed using SC networks, if the op-amp slew rate is not taken 
into account. 

The design approach used in this research, in general, is the extension of the 
design explained in [Ref. 6]. It addresses three different aspects of programmability. 
These are 

1. Programming the filter topology using a minimal set_of elements to obtain any 


f 
Dies of filtering function desired, for example LP, HP, BP, N and AP, by using 
a binary control word. 


tN 


Progamming the filter’s transfer function pganeicters (pole resonant frequency 
and a et Or ) for a chosen type of filtering function by using a digital 
control word. 


lsd 


pee eens Ue filter's transfer function (eens by changing the clock 
requency (practically from 50 KHz to 250 KHz). 
B. THE PROPOSED GIC SC PROGRAMMABLE FILTER 

The basic active network considered as the heart of the GIC SC programmable 
filter is the Generalized Immittance Converter (GIC) structure [Ref. 7] of Figure 3.1, 
whose superior performance was established in the [Ref. 8] and [Ref. 9]. The filters 


transfer functions derived can be found in [Ref. 6]. 





Figure 3.1 The GIC Implementation Using Op. Amps.. 


The bilinear realization of a continuous resistance was used throughout the 
design. Table 1, illustrates that for any filter realization, five switched capacitors, two 
unswitched capacitors, and two op. amps. are required. G stands for the admittance 
represented by switched capacitor. 

The passive elements are connected to the different nodes, shown in Figure 3.2, 
for the different realizations. A set of MOS bilateral switches controlled by a digital 
binary word are used to interchange the elements to achieve the different types of filter 
realization. The truth table of the switch control logics 1s shown in Table 2. Figure 3.3 
illustrates the CMOS logic circuit used for realizing this truth table. 

Four of the five switched capacitors are equal and of value Ca each, the fifth 
switched capacitor is the oF determining capacitor and of value C, =C ear Other 
two capacitors in the circuit are equal, and of value C = dt ail’, each. The two 
equal banks of capacitors are used to control Wo Each bank contains n binary 
weighted capacitors connected in parallel through analog CMOS switches as shown in 
Figure 3.4. Using a digital binary word of n bits to control Wo 2" different values of C 
will result at the 2 terminals of both capacitor banks that correspond to 2” different 


values of Wo: 
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Using a similar technique the value of C, can be controlled through a bank of m 
binary Weighted capacitors in parallel through analog CMOS switches as shown in 
Figure 3.5. Using a digital binary word of m bits to control Q. 2™ different values of 
C, can be achieved that correspond to 2™ different values of Q.: Thus, full 
independent control of the pole pair ve and Q, are achieved by programming the 


switches to obtain the corresponding C and C,. 


C. THE REALIZED GIC SC PROGRAMMABLE FILTER 

The detailed diagram of the signal processing circuit of the constructed GIC SC 
filter is shown in Figure 3.6. The value of m and n were selected m=n=4. Thus, 15 
different values of W, (f,) and Q, were obtained as it is illustrated at the corresponding 
Table 3 and Table 4. The implemented banks for the control of We and Qo are shown 
in Figure 3.7 and Figure 3.8 with their values. The clock frequency is 10! KHz 
throughout the design. According to this clock frequency, the switched capacitors value 
were choosen within range from 500 pF to 8nF, using experimental result of Chapter | 


Rigune lel>- 
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Figure 3.2 Schematic Diagram of the Programmable GIC SC Filter 
Showing the Controlled Nodes and Switches. 
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Figure 3.5 The Capacitor Bank for the Programming of Q,. 
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IV. COMPUTER SIMULATION OF 


THE PROPOSED GIC ACTIV 
FILTER hie 


A. INTRODUCTION 

In order to observe the theoretical frequency responses of the different 
realizations, that is LP, HP, BP, N and AP, with different pole frequencies and quality 
factors, the computer program, GICFRQ, has been written in FORTRAN 
programming language. The Display Integrated Software System and Plotting 
Language (DISSPLA) has been used to produce the desired graphic output. This 
frequency response program is shown in Appendix B. As it is seen this program, 
accepts some inputs from a data file, related with the transfer function of the 
realization. These inputs are processed through the program and results are passed to 
DISSPLA to draw the graph. One can plot the phase responses of the realizations as 
Well as the magnitude responses by making minor changes in the DISSPLA package. 

Because of the two four-bit control words for each realization, there will be a 
rather large number of input sets. For this reason a few intermediate control words 
were picked up randomly. In order to compare the computer simulations to the 


experimental results, the same control words are going to be choosen in next chapter. 


B. SIMULATION RESPONSE(S) 
1. Low Pass Filter Realization 
Using the elements values prescribed in Table 1 yields the following low pass 


filter transfer function 


2w? 


H(s) = 5 ) 5 
S BN Q,)s -'- en 


(eqn 4.1) 


First, it is assumed that the control word for Q, is set to 1010 which corresponds to Q, 
= 1.66 as shown in Table 4. Then Figure 4.1 through Figure 4.3 are plotted by varying 
the Wale ¥ to illustrate the magnitude responses of the ideal LPF GIC realization. 
Second, the control word We is set to 1010 which corresponds to i. = 2.595 KHz as 
shown in Table 3. Then Figure 4.4 are plotted for oF = 2.11, Q, = 4.07, Q, = 5.192. 
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Figure 4.2 Ideal GIC LPF Magnitude Response. 
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Figure 4.3 Ideal GIC LPF Magnitude Response. 
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Figure 4.4 Ideal GIC LPF Magnitude Response with Different QS: 
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2. High Pass Filter Realization 
Using the elements value shown in Table | yields the following high pass filter 


transfer function 


25° 


GO pert ac 
s7+(W/Q.)s+ W, 


(eqn 4.2) 


If we do the same set up as we did for the low pass filter case, we can obtain the 
various magnitude response curves shown in Figure 4.5 through Figure 4.7. Figure 4.8 
illustrates the effect of changing Q, on filter response. The same quality factor was 
especially chosen as in the case of the low pass filter, because similarities and 


differences can be illustrated more apperant in this fashion. 
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Figure 4.5 Ideal GIC HPF Magnitude Response. 
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Figure 4.6 Ideal GIC HPF Magnitude Response. 
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Band Pass Filter Realization 


U2 


The band pass transfer function of the GIC filter is 


2(W_/Q_)s 
tlle ceteeeaneeeeanEEEeee. on (eqn 4.3) 
s- + (W /Q,)s ~~ ee 


This will correspond to a gain of two at the pole frequency 

HW) = 2 
BPF realization is simulated for various Wade) when oO = 4.52 in Figure 4.9 and 
Figure 4.10. Figure 4.11 shows the magnitude response of BPF with different qualitv 
factors namely Q, = 5.192, Q, = 7.796 and Q, = 15.54 with Bs held constant. All 
figures show that the maximum gain at the pole frequency (W) is 2 independent of the 


quality factor. 
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Figure 4.9 Ideal GIC BPF Magnitude Response. 
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Figure 4.10 Ideal GIC BPF Magnitude Response. 
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Figure 4.11 Ideal GIC BPF Magnitude Response with Different QS: 


4. Notch Filter Realization 
Using apropriate element combination, a notch filter transfer function can be 


obtained as 


se+W 2 
n 


74 +W? 
7+ (WQ,)+W, 


H(s) = 


Figure 4.12 illustrates the ideal notch filter amplitude response for a variety of 
frequencies with constant Q, Oe 1.66, while Figure 4.13 for a variety of Q,$ with 
constant W (fo): The notch frequency, W_, is determined by both pole frequency and 
quality factor.As it is shown in the figures, if one of the parameters is held as constant 


the notch frequency still will continue to shift along the frequency axis. 
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Figure 4.12 [deal GIC NF Magnitude Response. 
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Figure 4.13 Ideal GIC NF Magnitude Response with Different Q,S: 
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5. All Pass Filter Realization 
As proposed in Table | using the same element values as in the notch filter 


case but different output node, all pass filter transfer function can be derived as 


s*-(W_'Q )s+W? 
ee 


AS 3 D 
2+(W_'Q )st+W? 
S or Q,)s We 


(eqn 4.35) 


which takes the value of Magnitude - | as ideal amplitude response. All pass transfer 
functions are often needed for delay equalization, that is, the phase response of the 
circuit has to have 360° shifting property around the pole frequency, . This property 


agree with the computer simulation result shown in Figure 4.14. 
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Figure 4.14 Ideal GIC AP Phase Response. 
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V. EXPERIMENTAL RESULT_OF THE PROGRAMMABLE GIC SC 
FILTER 
A. GENERAL 
All parts of the circuit, of which detailed designs were shown in Chapter III, were 
integrated to construct the complete digitally controlled programmable GIC SC filter. 
After the circuit was built a variety of measurements were taken in order to study the 
response of the network with the different configuration, as explained in Chapter III. 
The switches in this GIC SC filter, are devided into two sets, each performing a 
different function. The first set of switches are used to control the filter characteristics. 
The second set are used for resistor realization using the switched capacitor techniques 
discussed earlier. The effect of the control switches which introduce a resistance of 80 
ohm each at closed position can be found in [Ref. 6]. There are some general problems 
which may cause changes in the response of the network. These are listed below. 
le eee board creates capacitance values in the order of pfs at each 
2. The master clock was provided by the EXACT MODEL 120 waveform 
generator which 1s not so stable. Its output varies with the time. 


3. All switches introduce a resistance value from 60 ohm to 90 ohm at closed 
position. 


All probes being used for measurement have some capacitance value. 
Grounding problem. 


The equipment scaling problem exists during the measurements. 


eee 


Capacitor values change within their tolerance limitation. 
In order to compare the experimental result to the computer simulation graphs, 
the closest control bitwords were used to set up the pole frequency and quality factor. 
After the response of all filter configurations were plotted, the clock frequency was 
varied to observe changes of pole frequency and the amount of linearity between clock 
and pole frequency. Finally, measurements were taken to show the sensitivity the 


quality factor to clock frequency changes. 


B. LOW PASS FILTER 
With the topology-control bitword 000, the network realized a Low Pass Filter 
response. Altough all calculations were made under 101 KHz clock frequency, because 


of the unstability of the waveform generator mentioned above, average clock frequency 
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was found to be 99.5 KHz. Q, was set to 1010 which corresponds to approximately 
1.66. Magnitude response curves were plotted within the range of programmability of 
pole frequency, namely from 1.81 KHz to 30.3 KHz. Expected results were obtained 
up to 7.55 KHz. Beyond this frequency, the gain declined slightly, while quality factor 
experienced little changes. This can be seen in Figure 5.1 through 5.3. Control bitword 
for pole frequency was set to 1010, which is f,=2 KHz, while changing Q, control 
bitword to 1000, 0101 and, 0011 which corresponds to 2.11, 4.07 and, 5.19 respectively. 
Figure 5.4 was plotted under these conditions. The plot is quite similar to computer 
simulation except Q, = 9-19. This is also expected because of the fact that setting up 
the quality factor to 5.19 requires 440 pf capacitance value. This value falls out of the 
range that we considered in Figure 1.15. The continuous resistance value was not 
realizable according to the design value. As a result, some deviation from quality 
factor was experienced. Figure 5.4 shows that the pole frequency was not effected, 


since the 440 pf capacitor is not the element used to control to 
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Figure 5.1 Experimental Result of LPF with Different tS 
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Figure 5.3 Experimental Result of LPF with Different f 5 
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Figure 5.4 Experimental Result of LPF with Different Q,. 


C. HIGH PASS FILTER 

With the topology-control bitword 001, the network realized a High Pass Filter. 
For this case, all measurements are taken under 101.8 KHz clock frequency. The 
qualitv factor was set to 1.66 for six different values of pole frequency shown in Figure 
5.5 through Figure 5.7. In regard to pole frequency and gain, expected results were 
obtained. But the quality factor tends to get closer to computed curve up to 5 KHz 
frequency range and then diverges. On the other hand, if we consider changes in the 
quality factor, as shown in Figure 5.8, the small Q,3 were implemented more accuretly. 
This 1s because, small quality factors require larger capacitance value and these values 
can realize the equivalent resistance Well enough. Three different quality factors were 
used when Figure 5.8 was plotted. There is no distortion in the pole frequency as 


expected theoretically. 
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Figure 5.6 Experimental Result of HPF with Different f + 
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D. BAND PASS FILTER 

The 010 topology-control bitword realizes the Band Pass Filter. The quality 
factor was set to 0100 which corresponds to 4.52. Figure 5.9 and Figure 5.10 were 
plotted for various pole frequencies. The results agree with the computer simulations 
with minor amplitude fluctations. Figure 5.11 was poltted at 2.59 KHz pole frequency 
for various quality factors, 5.19, 7.76, 15.5. An important observation on this figure is 
the fact that the quality factor increases, while the gain decreases. But again, similar to 
other cases, pole frequency remained constant. This situation is shown apperantly in 


Figure Sage 
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Figure 5.9 Experimental Result of BPF with Different 
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Figure 5.11 Experimental Result of BPF with Different Q,: 
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E. NOTCH FILTER 

With the topology-control bitword O11, a Notch Filter realization can be 
achieved. Figure 5.12 illustrates the amplitude response for different frequencies for 
Q, = 9-19, while Figure 5.13 illustrates amplitude response for fa 755 KHz and a 
variety of quality factors. The above figures agree with the theoretical transfer 
function. Since the notch frequency, W_, depends on the both pole frequency and 
quality factor, whenever pole frequency (or quality factor) changes, notch frequency 
and also quality factor (or pole frequency) reflect this variation as it is seen in Figure 


5.12 and Pigiitievos: 
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Figure 5.12 Experimental Result of NF with Different fe 
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Figure 5.13. Experimental Result of NF with Different Q, 


F. ALL PASS FILTER 

The topology-control bitword 100 realizes an All Pass Filter. Figure 5.14 
illustrates frequency response of the realized all pass filter for different pole frequency 
while quality factor is held constant. For this particular filter, 360° phase shift was 


observed at the pole frequency which agrees with the computer simulation. 


G. EFFICIENCY OF THE CLOCK FREQUENCY 

In order to observe the affect of the clock frequency, f,, two different plots have 
been taken. As mentioned in Chapter III, if the clock frequency changes by some 
constant @, the pole frequency shifts by an amount which is related to this constant @. 
The programmable filter can be programmed to almost any frequency (from | KHz to 
60 KHz) by using this feature. Figure 5.15 shows this property, where the first LPF 
uses 50.97 KHz while the second is using 103.5 KHz. This means that the clock 


frequency changed almost 50 KHz, correspondingly pole frequency shifted almost 0.75 


IE 


KHz. The last three LPF curves in Figure 5.15 had 103.5 KHz, 201 KHz and, 299.5 
KHz clock frequency. The difference from each other was almost 100 KHz, 
correspondingly the pole frequencies differ from each other almost 1.5 KHz which 
proves the concept stated above. 

Finally, the affect of the clock frequency with respect to quality factor, Q, will 
be demonsrated. In order to do so, first pole frequency was set to 1111, 1.81 KHz by 
using 103.25 KHz clock frequency (curve 1 in Figure 5.16) then pole frequency, to was 
changed to 4.25 KHz using control bitword 0111 with same clock frequency (curve 2 in 
Figuré S216), From this to the clock frequency was decreased until the same ie was 
realized as in curve 1. This was plotted as curve 3 in Figure 5.16. By comparing curve 1 
and curve 3 one can observe that there is no major changes in terms of quality factor. 
This can prove experimentally the low sensitivity of the quality factor Qo lo the cide. 


frequency variations. 
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Figure 5.14 Experimental Result of AF with Different fo 
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VI. CONCLUSION 


The design described here has resulted in a universal programmable switched 
capacitor filter that can be digitally controlled to realize almost any practical filter 
specifications. This is done through the use of CMOS switches controlled by binary 
codes to program the filter topology, the filter center frequency, quality factor. The 
bilinear switched capacitor realization of a continuous resistor was used as a key 
concept of the research. The design procedure required developing optimum switching 
arrangments for the minimum redundancy in components and least dependence of the 
filtering function on switching imperfections such as switches stray capacitances and 
non-zero and nonlinear switch-on resistance. Further investigation 1s needed eliminate 
the effect of the stray capacitance introduced by the switched capacitors. The 
sensitivity of 0; is found to be low with respect to clock frequency. The clock 
frequency could also be used to program the filter by scalling the center frequencies 
through different ranges. The experimental results showed close agreement between 
theory and practice. Further, these results indicate that these realizations are insensitive 
to temprature and power supply variations. 

This study can be extended for developing a wide bandwidth programmable 
switched capacitor filter using the composite operational amplifier thechnique proposed 
by [Ref. 9]. Such implementetation would lead to a very useful monolithic device at 
moderate cost. 

This research can be also brought into the subject of VLSI design to implement 
the whole network structure into the single custom chip, which is the final goal of the 


research. 
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APPENDIX A 
TRANSFORMATION PROGRAM 


SSTORAGE:Z 


SPAGE 


101 
OZ 
103 
104 
105 
106 
107 
108 


RESew ool Ac, 60,6L,82,K,F,A,SP,DELTH 
REAL AZO,AZ1,AZ2,BZ0,BZ1,B22 


SP = 0.0 
DELTH = 0.0062831853 
Me 2 
= 
NP = 501 
BZ0 = 1.0 
WRITE(*,*)'THIS PROGRAM TRANSFORMS SECOND-ORDER S-DOMAIN TRANSFER’! 
WRITE(*,*)'FUNCTION TO Z-DOMAIN TRANSFER FUNCTION USING BILINEAR' 
WRITE(*,*)'MAPPING METHOD.WT/2 << 1 IS ASSUMED.' 
WRITE(*,*)'IF YOU ARE DEALING WITH THE SC CIRCUIT,YOU HAVE TO! 
WRITE *,*)'TAKE SAMPLING FREQUENCY AS TWICE AS MUCH ITSELF. ' 
WRITE(*,*)'ENTER THE SAMPLING FREQUENCY IN HZ.'! 
READ ( ne F 
WRITE(*, i. (AO*S**2)+(A1*S)+A2 IS THE FORM OF NUMARETOR' 
WRITE(* x ENTER AO,Al,A2' 
READ ( Ne AO,Al1,A2 
WRITE(*, f (BO*S**2)+(B1*S)+B2 IS THE FORM OF THE DENOMINATOR' 
WRITE(*,*)'ENTER BO,B1,B2' 
READ(*,*) BO,B1,B2 

= 2.0%F 

= (B2+ (BL*K) + (BO*K*A2 ) 
AZO = (A2+(Al ie GE AG 2))/A 
RZie= (2. 9*A2) -(2.0*A0*K /B 
AZ2 = (A2-(AL*RK)+(AO*XK*AZ 
BZ1 = ((2.0%B2)-(2.0%(BO*K KD)))/R 
BZ2 = (B2-(B1*K)+(BO*K**2))/A 
OPEN(8,FILE='BILINEAR.DAT' , STATUS='NEW' ) 
WRITE(8,101) N,M,NP,SP,DELTH 
WRITE(8,102) AZO 
WRITE(8,102) AZ1 
WRITE(8,102) AZ2 
WRITE(8,102) BZO 
WRITE(8,102) BZl 
WRITE(8,102) BZ2 
CLOSE (83 
WRITE(*,*)'H(Z) = ((AZO*Z**2)+(AZL*Z)+AZ2) / ((Z**2)+(BZ1*Z)+BZ2) ! 
WRITE(*,101) N,M,NP,SP,DELTH 
WRITE(*,103) AZO 
WRITE(*,104) AZ1 
WRITE(*,105) AZ2 
WRITE(*,106) BZO 
WRITE(*,107) BZ1 
WRITE(*,108) B22 
FORMAT (213, 2X,14,2(4X,F13.10)) 
FORMAT (3X,F13.103 
FORMAT(3X,'AZO = ',F13.10 
FORMAT(3X,'AZ1 = ',F13.10 
FORMAT(3X,'AZ2 = ',F13.10 
FORMAT (3X,'BZO = ',F13.10 
FORMAT(3X,'BZ1 = ',F13.10 
FORMAT (3X,'BZ2 = ',F13.10 
STOP 
END 
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APPENDIX B 
SIMULATION PROGRAM 


THIS PROGRAM CALCULATES THE FREQUENCY RESPONSE OF A CONTINUOUS ... 
. SYSTEM. THE INPUT IS A RATIO GP THO Perici rssOr inh es OnTit.. 


(B(O)*S*4M + B(1)*S**(M-1) + ... + B(M-1)*S + B(M)) / 
(A(O0)*S*AN + A(1)*S**(N-1) + 2... + A(N-1)*S + A(N)) 


WHERE THE NUMERATOR COEFFICIENTS B(0), Be oa BM) ANDETHE 


THE DENOMINATOR COEFFICIENTS A(O), A(1 ),.0+,A(N) ARE REAL ... 
NUMBERS AND THE DEGREE OF THE NUMERATOR POLYNOMIAL, M, Ds 
THE DEGREE OF THE DENOMINATOR POLYNOMIAL, N, ARE POSITIVE ° 
INTEGERS. THE A'S,B'S,M, AND N ARE READ IN AS DATA, AS IS cee 
OMEGAO (THE STARTING VALUE OF OMEGA IN RADIANS rene 
OMEGA IS DEFINED BY S = J*OMEGA DLOMGA (THE 

INCREMENT OF OMEGA), AND NUMPTS (THE NUMBER OF POINTS) 


IT IS ASSUMED THAT ALL SYSTEM POLES HAVE NEGATIVE REAL PARTS. 


INTEGER M,N,NUMPTS 
COMPLEX S,DEN,NUM,H,CI 
REAL MH(1000) ,PH(1000) , OMEGAV(1000),A(128),B(128) , IMS,RES, DLOMGA 
CALL TEK618 
CALL SHERPA (° TKISERI ' tat 3) 
CALL BLOWUP(0 75) 
CALL PAGE(9.0, aaa) 
CALL HWROT('AUTO! 
CALL NOBRDR 
CALL HEIGHT(0.15) 
CALL MX1ALF('STANDARD' , '&') 
CALL MX2ALF 'L/CSTD', '#!) 
CALL MX3ALF('GREEK' ,'%') 
CALL MX4ALF('L/CGREEK', '@') 
NP = 0 
NPROP = 1 
READ(4,1000) M,N,DLOMGA, OMEGAO ,NUMPTS 
NP = NP+1 


WRITE(8,1008 
WRITE(8,1110) M 
WRITE(8,1111) N 
WRITE CS ZU Mi nS 
WRITE(8,1113) OMEGAO 
WRITE(8,1114) DLOMGA 
IT IS ASSUMED THAT DLOMGA IS READ IN IN RADIANS/SEC 
PI=3.141592654 


MP1=M+1 

NP1=N+1 

See elo Bae 
READ(4,1001 A(I),1=1,NP1 
WRITE(8,1004 

WRITE(8,1005) (B(I),I=1,MP1) 
WRITE(8,1006 

a a “HO (A(I) ,I=1,NP1) 


BEGIN MAIN LOOP 


DO 100 L=1,NUMPTS 
NUM = CI*B(1) 
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DEN = CI*A(1) 
E UPDATE VALUE OF OMEGA 
OMEGA = OMEGAO + {usd *XDLOMGA 


OMEGAV (L) = OMEGA/(2.*PI) 
RES = 
IMS = OMEGA 


S = CMPLX(RES,IMS) 
IF(M. ae or COG wae 
DO 50 vil 
50 NUM = s*NUM re CI*B(J+1) 
60 CONTINUE 
IF(N. ee .0) GO TO 80 
DO 70 = 1,N 
70 DEN = S*DEN + CI*A(J+1) 
80 CONTINUE 
H = NUM/DEN 


MH(L) = CABS(H 
IF oy sGEaie ae GO - 90 
IF(ABS(IMAG(H)).LE.1.E-15) PH(L) = 0.0 
IF (IMAG a ‘cr. 7 neo Ee = ao 
IF(IMAG(H).LT.-1.E-15) PH(L) = -PI/2. 
Comino ¢ 

90 PH(L) = ATAN2 ee ,REAL(H) ) 

95 PH(L) = PH(L)*180 


100 CONTINUE 
WRITE( G7 LCC9 
WRITE(8,1002 
Preise (oO 7 
ae) = i NUMETS 
150 WRITE(8, 1003 )OMEGAV(L) , eae Pench) 
1000 FORMAT ‘ ya Z(ElGe4) 13 
El 


1001 FORMAT 
1002 FORMAT / , 8%, JOMEGA?, nN, Eee Eeciay i ae jon, PEASE CP.) 


2 
: 
1007 FORMAT(8X,'RAD/S ',23X,'DEGREES' 
1003 FORMAT(3(2X,E14. 6))’ 
1004 FORMAT(2X,'THE NUMERATOR COEFFICIENTS B(0),B(1),...,B(M) ARE',/) 
1005 FORMAT(10(1 x, mi:4).//) 
1006 FORMAT(//,2X, 'THE DENOMINATOR COEFFICIENTS A(0),A(1),...A(N) ARE', 
1008 FORMAT(20X,' INPUT DATA ',///) 
1009 FORMAT(//////,20X,' OUTPUT DATA ',///) 
1010 FORMAT(2(2X,E11.4,/ 
1110 FORMAT(2X,'DEGREE OF NUMERATOR = ',I3) 
1111 FORMAT(2X,'DEGREE OF DENOMINATOR = ',13) 
1112 FORMAT(2X,'NUMBER OF FREQUENCY POINTS = er 
1113 FORMAT(2X,'STARTING VALUE OF OMEGA = ',E12.6 
1114 FORMAT(2X,'INCREMENT OF OMEGA = ',E12.6,/) 


KAKKKAKKKKKK DISSPLAY *AAAKAAAAAKKK 


CALL COMPRS 
XMIN = OMEGA 
XMAX = OMEGAY (NUMPTS) 
XMAX = 5000. 
XINC = XMAX/20. 
MAGMAX = -0.05 
MAGMIN = 0.0 
PHMAX = hee 
PHMIN = 1.E15 
DO 180 fe = 1, NUMPT 
IF(MH(1I).GT. MAGMAX ) a aeuan = MH(I) 
oe PH(I).GT.PHMAX) PHMAX = PH(I) 
F(MH(I).LT. PN), MAGMIN = MH(I) 
180 TF(PHOLY. Ee PHMIN) P HMIN = PH(I) 
GO TO 2 10, 220,220,220) ,NP 
GO Tare lon 72) 220,220,220 ,220),NP 
210 CONTINUE 
CALL AREA2D (6. 03.5) 
GALL ae FREQUENCY H#Z&S'! 
CALL YNAME('MAGNITUDE $', 1003 


Qa ANN 


100) 
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220 


CALL HEADIN(' $',100,1.2,4 
CALL HEADIN(' $',100,1.2,4 
CALL HEADIN(' Rl = 21.77K C2 = 431#P&F R3 = 1.17M$',100,1.2,4) 
CALL HEADIN(' FREQUENCY RESPONSE MAGNITUDES! ,100,1.2,4) 
CALL CROSS 
CALL XNONUM 
CALL GRAF(XMIN, 'SCALE' ,XMAX,0.,'SCALE' ,MAGMAX) 
CALL GRAF (KMIN, "SCALE! ,XMAX,0.,0.5,2.5) 
X1=MAGMAX/14 
CALL HEIGHT (0.2 
CALL RLMESS S58 | eloe 3290.0) 
CALL HEIGHT 5) 
CALL RLMESS ee ote 3 
CALL RLMESS('-3@P&/45',100,-2.60,-X1 
CALL RLMESS('-@P&/2$', 200 1214 
CALL RLMESS('-@P&/4$' ,100,-0.98,-X1 
CALL RLMESS('0$',100,-0.037,-Xl 
CALL RLMESS('@P&/43S', en 2 -6o aa 
CALL RLMESS('@P&/2$! ,100,1.47,-X1 
CALL RLMESS('3@P&/4$' ,100,2.22,-X1) 
CALL RLMESS('@P&$! ,100,3.10,-X1) 
CALL MARKER(15) 
CALL CURVE (OMEGAV,MH,NUMPTS ,0) 
CALL AREA2D(6.0,3.5) 
CALL KNAME ("FREQUENCY Hi HEZES ', 100) 
CALL YNAME('P#HASE&$',100) 
CALL HEADIN('  $', ot: eo 
CALL HEADIN(' §$',100,1.2,4 
CALL HEADIN(' LOSSLESS INTEGRATOR $',100,1.2,4) 
CALL HEADIN('R1 = 21.77K C2 = 431#P&F RS = L.17H$' 100, ea) 
CALL HEADIN(' Glc AP ECOu—nincey POOn tsa 
CALL HEADIN(' FREQUENCY RESPONSE PHASE asi 100,1.2) 
CALL CROSS 
CALL GRAF(XMIN, XINC, XMAX, PHMIN, 'SCALE! , PHMAX) 
CALL GRAF(XMIN,'SCALE! , XMAX, PHMIN, 'SCALE' , PHMAX) 
X2=(PHMAX- “PHMIN Vick 
CALL HEIGHT(O. 
CALL RLMESS 1305! ,100,3.2,0.0) 
CALL HEIGHT(0.15 
CALL RLMESS('-@P | 100 -saadaane 
CALL RLMESS('-3@P&/4$' ,100,-2.60, -X2 
CALL RLMESS ms ee 
CALL RLMESS(' 
CALL RLMESS('0$' ,100,-0.037,-X2 
CALL RLMESS('@P&/4$' ,100, aie 
CALL RLMESS( '@P&/25!' ,100,1.47,-X2 
CALL RLMESS '3@P&/4$! , 100, 2.22, -X2) 
CALL RLMESS('@P&$',100,3.10, -X2) 
CALL CURVE (OMEGAV, PH,NUMPTS,0) 
CALL RESET('CROSS' ) 
CALL ENDPL(0) 
IF(NP.LT.NPROP) GOTO 1 
CALL LINES('1.81KH#Z&$' , IPAK, 1 
CALL LINES('1.93KH#Z&5' , IPAK, 2 
CALL LINES('2.87KH#Z&5' , IPAK, 3 
CALL LINES('3.17KH#Z&5' , IPAK, 4 
XW = Beeb TEA as 
YW = YLEGND(IPAK,4 
XL = 6.2-.5- “KH. 
YL = 2.3-.5-YW 
@rGist SiNKI(aLe. 12 ner 2-.5,YL-.1,2.3-.5,2) 
CALL DOT 
CAmmecernCiee 1) 
ean a 
CALL RESET('BLNK1') 
CALL LEGEND(IPAK,4,XL, YL) 
CALL ENDPL(O) 
CALL DONEPL 


84 


AANAANANANANAAAANAAANANAANAAANANAANANAAANANNAAANANAAANANANNAANAANAANAAAAAANANAAANAANAANAANAANANANANANAANANANN 


neo 
Dae 


KAEAEKKKKKAKKKKKKKKKAKKKKKKKRKAKKRKRKKAKKKKAKKKKAKKK 


s INSTRUCTIONS FOR USE OF PROGRAM . 
KAKKKKKKAKKKKKKKKKKKAKAKKAKKAKAKKKKKKAKAAAK 


SER'S INPUT DATA IS TO BE LOCATED IN A FILE WHOSE FILENAME FILE- 
AND FILEMODE ARE FILE FTO4FO0O1 Al , THE OUTPUT DATA WILL BE 


STORED IN THE FILE FILE FTO8FOO1 Al. 
Se pe ee Sooo ree vee sooo ores oes SS RSS e595 
$ TO COMPILE, LOA ECUTE THE PROGRAM ONCE THE INPUT DATA FILE 
$ HAS BEEN CREATED, “TYPE <FROM FLIST> 
RU 
: THIS IN INVOKES THE FORTRAN VS COMPILER, LOADS AND EXECUTES THE 
$g5$55S5995$SS$$$$$$S$$5$9$$$S$$$S$$5S$$$S5559$$59$55SS$59S$S9S5S855 
THE INPUT DATA REQUIRED ARE (NOTE:REAL NUMBERS WITH DECIMAL POINTS. 
MUST BE IN THE FORMAT X.Y, EVEN IF X,OR Y, OR BOTH ARE ZERO): 
NAME TYPE RANGE(IF ARRAY) RESTRICTIONS (TF ANY) 
M INTEGER QO <= 4 12 
N INTEGER 0 <= N = 133 
OMEGAO REAL 
DLOMGA REAL 
NUMPTS REAL 1 <= NUMPTS <= 1000 
B(L} REAL 0,1, ...,M 0 <= M <= 128 
A(L REAL O,1, ...,N 0 <= N <= 128 
WHERE ... 
M = DEGREE OF NUMERATOR POLYNOMIAL 
N = DEGREE OF DENOMINATOR POLYNOMIAL 
OMEGAO = STARTING VALUE OF OMEGA (AS IN S = J*OMEGA)) IN RAD/S 
DLOMGA = INCREMENT OF OMEGA IN RADIANS 
B(L) = cake Ge IN NUMERATOR POLYNOMIAL (IN ORDER B(0),B(1), 
Bae) = oe a ai IN DENOMINATOR POLYNOMIAL (IN ORDER A(0),A(1), 
11. ACM 
THE FOLLOWING INPUT DATA SET PROVIDES THE FREQUENCY RESPONSE FOR THE 
SYSTEM CHARACTERIZED BY 
H(S) = 10S/(S+1)(S+5)/OMEGA 
FOR OMEGA = 0 TO OMEGA = 4.0 IN STEPS OF OMEGA = 0.2 
001 002 OZ 0.0 021 
10. 0.0 0.0 
1.0 6.0 5.0 
INPUT DATA 
DEGREE OF NUMERATOR = 1 
DEGREE OF DENOMINATOR = 2 
NUMBER OF FREQUENCY POINTS = 21 
STARTING VALUE OF OMEGA = 0.000000E+00 
INCREMENT OF OMEGA = 0.200000E+00 
THE NUMERATOR COEFFICIENTS B(0),B(1),...,B(M) ARE 
0.1000E+02 0.0000E+00 
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THE DENOMINATOR COEFFICIENTS A(0),A(1),...A(N) ARE 
0.1000E+01 


OVDDDADAWDADAWMDAIVWDIAIAIOMDAMAIOO0000 


OMEGA 
RAD/S 


.0OO0000E+00 
-200000E+00 
-400000E+00 
-600000E+00 
- 800000E+00 
- LOCCOGCE+TO!L 
-L120000E+01 
. L40000E+01 
. LEO000E+01 
- L8QC00E+01 
. 2O0O000E+01 
.220000E+01 
-240000E+01 
-2600000E+01 
-280000E+01 
. 300000E+01 
- 32Z0000E+01 
. 340000E+01 
. 360000E+01 
. 380000E+01 
-400000E+01 
STOP 

END 


0.6000E+01 


0.5000E+01 


OUTPUT DATA 


MAGNITUDE (M) 


.O0O0000E+00 
.391919E+00 
-740416E+00 
.102166E+01 
P1233 70 FCT 
-138675E+01 
-149402E+01 
-156719E+01 
.161531E+01 
.164497E+01 
.166091E+01 
.166654E+01 
.166435E+01 
.165616E+01 
.164335E+01 


tie 
OVDWDAIDODVAIO9O000000 


DOVDVDOVDVDOVWMDAOAIODOVAOOOVO900000O 


> See Onl =@ 
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PHASE(P) 
DEGREES 


-OO0000E+00 
.763994E+02 
.636246E+02 
»521935E+02 
~-422499E+02 
.336901E+02 
.263098E+02 
.198954E+02 
.142607E+02 
.925574E+01 
.476366E+01 
.694467E+00 
~S0Z1LI3Et01 
.643690E+01 
.959499E+01 
.162698E+01 -0. 
.160786E+01 -Q. 
566658701 -O. 
.156386E+01 -0. 
~153990E+01 -0. 
.246236E+02 


125288E+02 
152652E+02 
178262E+02 
202298E+02 
224913E+02 


Ul 
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